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Abstract 

The Born-Oppenheimer potential for the X S+ state of H2 is obtained in the range of 0.1 - 20 au, using 

(N ; 

analytic formulas and recursion relations for two-center two-electron integrals with exponential functions. 
For small distances James-Coolidge basis is used, while for large distances the Heitler-London functions 
Lj ■ with arbitrary polynomial in electron variables. In the whole range of internuclear distance about 1CP 15 

precision is achieved, as an example at the equilibrium distance r = 1.4011 au the Born-Oppenheimer 
potential amounts to —1.174475 931400 216 7(3). Results for the exchange energy verify the formula of 



PACS numbers: 31.15.ac, 31.50.Bc 



Herring and Flicker [Phys. Rev. 134, A362 (1964)] for the large internuclear distance asymptotics. The 
presented analytic approach to Slater integrals opens a window for the high precision calculations in an 
arbitrary diatomic molecule. 
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I. INTRODUCTION 

In order to accurately obtain vibrational and rotational spectra of molecules not only the non- 
relativistic Born-Oppenheimer (BO) potential has to be calculated, but also the adiabatic, nona- 
diabatic, relativistic and quantum electrodynamics effects. An excellent agreement of theoretical 
results yj] with experimental values achieved for dissociation energies of the hydrogen molecule 
H 2 D20 and D 2 [3D indicates good understanding of all physical effects up to 1CT 8 precision level. 
Such a precision has not yet been achieved for molecular systems with more than two electrons. 
It is because numerical calculations with explicitly correlated functions are very demanding. The 
commonly used explicitly correlated Gaussian functions require global minimization of thousands 
of nonlinear parameters, and the accuracy achieved for relativistic effects is difficult to estimate, 
due to the improper analytic properties of these functions. The largest system considered so far is 
the helium dimer, where about 1CT 5 precision is achieved for the interaction energy Q4£] . 

Herewith, we apply explicitly correlated exponential functions (for a recent review see Ref. 
|5|]), in order to increase numerical accuracy obtained so far for two electrons diatomic systems. 
The use of exponential functions may allow for a significant improvement not only in rotational 
and vibrational energies, but also in other important properties like magnetic shielding, spin- 
rotational or spin-spin coupling constants. It is because evaluation of relativistic effects is very 
sensitive to correct analytic properties of the basis functions, a good example being the relativistic 
correction to the magnetic shielding |g). 

When 1CT 10 precision is achieved for vibrational transitions, the electron-proton mass ratio 
can be determined with better accuracy, than it is known presently [7|]. Moreover, any non- 
electromagnetic long-range interactions between nuclei, which cannot in principle be excluded, 
may be visible not only in vibrational spectra, but also in the spin interaction between nuclei [8Q, 
for example in HD molecule. The electromagnetic spin interaction is extremely small, of order 
m 3 /(m p rrid) a 6 , about 43 Hz for the ground molecular state [9], and any deviation between exper- 
imental values and theoretical predictions will signal existence of a nonelectromagnetic long-range 
interactions between nuclei. 

In this work we present an effective computational method with the use of exponential functions 
and numerical results for the Born-Oppenheimer energy of the H 2 molecule. We demonstrate that 
1CT 15 numerical precision can be achieved for BO energy with about 22 000 basis functions using 
the analytic formulae and recursion relations for Slater integrals, which have been obtained in 



our previous work in II 1 QTI . The obtained BO potential is at least two orders of magnitude more 
accurate, than anything published so far, and at the equilibrium distance our calculations confirm 
the most accurate so far result of Cencek and Szalewicz in Ref. Ill ill well within their uncertainties. 
Our results for large internuclear distances verify the asymptotic formula of Herring and Flicker 



IJ12I1 . although we observe significant contributions from the off leading terms. Finally, we analyze 



possibility of the extension of this method to more than two-electron systems. 

II. SHORT RANGE OF INTERNUCLEAR DISTANCES 

The Schrodinger equation for H 2 molecule in the Born-Oppenheimer approximation is 

H(f)(fi,f 2 ) = E{r)4>{f u r 2 ), (1) 

where 

H = --i - -£ + — + -. (2) 

2 2 r 1A r lB r 2A r 2B r 12 r 

The most efficient basis set to solve this equation is the one introduced by Kotos and Wolniewicz 



in rtl3L_However at small nuclear distances r = r AB , we can use much simpler James-Coolidge 



basis LI 4], namely the functions of the form 



^c w (l + P4 B )(l + P 12 )e- Q 



(fu+ns) e -oc(r 2 A+r2B) r n 



12 



such that 



(riA ~ r 1B ) n2 (r 2A - r 2B ) n * (r lA + r 1B ) n ± (r 2A + r 2B ) n * (3) 



^rn<n (4) 



8=1 
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with Q = 3, ... 20. This basis have recently been used by Sims and Hagstrom 111511 for 10 
precision calculation of BO potential in the range r = 0.4 — 6 au, by a traditional way of numerical 
evaluation of the corresponding two-center integrals. Here, using Ref. [10] we derive analytic 
expression for all these integrals up to fi = 20. They have very simple form and involve only 
exponential integral (Ei), exponential and logarithmic functions. In this derivation we make use 
of a differential equation, which is satisfied by the most general two-center two-electron Slater 



integral, see Ref. fl 1 QTI 



For the presentation of analytic formulas we consider an integral 
f(ni,n 2 ,n 3 ,n4,n 5 ;r,u,w) = 



d 3 ri f d 3 ro e~ UTlA e~ ur ' 1B e~ WT2A e~ vr - l: 

47T 



4tt r XA r 1B r 2A r 2B r\ 2 ni 
{via - r 1B r (r 2A - r 2B ) n * {r lA + r 1B )" 4 {r 2A + r 2B )" 8 (5) 



at the internuclear distance r = 1, since arbitrary distances can be obtained from the relation 



/(m, n 2 , n 3 , n 4 , n 5 ; r, u, w) = /(m, n 2 , n 3 , n 4 , n 5 , l,ru,rw)r 



2+ni+n2+ri3+n4+n5 



(6) 



Positive powers of n 4 and n 5 can be obtained by differentiation with respect to u and w corre- 
spondingly, so without loss of generality one can assume n 4 = n 5 = 0. Let us introduce auxiliary 
functions yi defined by 



y = e~^ 

Vl = Ei(-2w) e w ~ u + Ei(-2u) e 

y 2 = Ei(-2w)e w - u -Ei(-2u)e 

y 3 = Ei(-2(u + w))e u+w + 



_ , 2uw 
ln[— — )+ 7 

u + w 

. . 2uw 

m[ — — j+7 

u + w 



-(u+w) 



-(u+w) 



(7) 

(8) 

(9) 

(10) 

(11) 



?/4 = Ei(-2 (u + w)) e u+w 

Then, using the equation (69) from IllOn . all nonvanishing integrals f{n\,n 2 ,n 3 ) = 
f(rii,n 2 ,n 3 , 0, 0; 1, u, w) with J2i n i < 4 are the following: 



/(0,0,0) 


/(I, 0,0) 


/(2,0,0) 


/(0,2,0) 


/(o,i,i) 


/(3,0,0) 


/(I, 2,0) 



Vz-Vi 
4\uw 
2/o 



Auw 



y Q (2uw + u + w) +y 2 (u-w) -y A (u + w) (y 3 - y x ) (u 2 + w 2 ) 



2Au 2 w 



2 n,\2 



2Au 3 w 



3 niyO 



Y2uw 



y y 1 (uw- 1) -y 2 (u -w)+y 3 (uw + 1) - y 4 (u + w) 



6uw ' 12 u 2 w 2 

y (u 2 w 2 + 3uw(u + w) + 3 (u 2 + w 2 )) 



24 u 3 w 3 



\2uw 



(12) 
(13) 

(14) 

(15) 

(16) 

(17) 
(18) 



y (2u 3 w 3 + 6u 3 w 2 + 9u 3 w + 3u 3 + Qu 2 w 3 + 4u 2 w 2 + 6m 2 w + 9uw 3 + 6-uw 2 + 3u> 3 ) 

J V ' ' / on. .4 4 



80m 4 u> 4 
1/2 (w — iu) (3w 2 + «w + 3u> 2 ) 1/4 (u + w) (3m 2 — «w + 3u> 2 ) 



40w 4 w 4 40w 4 u> 4 



(s/3 ~ S/i ) (« 4 w 2 + 3m 4 + w 2 w 4 + 2w 2 w 2 + 3w 4 ) y 3 + Vi (m 

40u 5 w 5 20u 2 w 2 

1/2 (w~ w)(uw + 2) _ y 4 (w + w)(ww-2) |/o(4mw + 3m + 3w + 4) 
/l ' ' j ~ 120w 3 u> 3 120w 3 u> 3 120w 2 u> 2 

yi (-u 2 + 2mu + w 2 — 2) y 3 (w 2 — 2-uw + w 2 — 2) 



/(2,1,1) 



120w 3 w 3 120w 3 u> 3 

y (2u 2 w 2 + 3-u 2 w; + 6m 2 + 3uw 2 + 6w 2 ) 

120u 3 w 3 
y 2 (u — w) [u 2 w 2 + 3-u 2 + 3u> 2 ) y A (u + w) (u 2 w 2 + 3w 2 + 3w 2 ) 



120w 4 w 4 120«% 4 

(ys + Vi) (u 2 + w 2 ) (y 3 - yi) (2u 2 w 2 + 3m 2 + 3w 2 ) 
40w 3 u> 3 120m 4 w 4 



(21) 



/(0,4,0) = %-^ (22) 

20 uw 



yijuw- 1) _ y 2 (u-w) y 3 (uw + 1) _ y 4 (u + w) _ y 
20u 2 w 2 20u 2 w 2 20u 2 w 2 20u 2 w 2 ~ lOuw 

y 2 (u — w)(uw — 3) y&(u + w)(uw + 3) y (uw + 2u + 2w + 6) 
15u 3 w 3 15u 3 w 3 15u 2 w 2 

2/3 + 2/1 (Vs ~ Vi) (3m 2 w 2 + 4u 2 + 4w 2 + 12) 
2u 2 w 2 60u 3 w 3 



Integrals with higher powers of m are similar, with simple polynomial form in 1/u and 1/w. We 
have tabulated all integrals, such that ^ n^ < 45 what is sufficient for matrix elements involving 
exponential functions with maximum value of CI = 20. 

III. LONG RANGE OF INTERNUCLEAR DISTANCES 

For the large internuclear distance the James-Coolidge basis is not appropriate as it does not 
include the Heitler-London function. Instead, we employ generalized Heitler-London functions, 
which are the product of a Heitler London function with an arbitrary polynomial in all electron 
distances, 

4> = J2 c {n} (l + PAB)(l + Pi2)e^ A+r ^r^r n 1 Xr n 1 : hrT A r^ (25) 

M 



such that 

5 

J> 4 <fi, (26) 

and call them the explicitly correlated asymptotic (EC A) basis. Matrix elements of the nonrela- 

tivistic Hamiltonian can be expressed in terms of direct integrals of the form 

[ d 3 n f d 3 r 2 r e~ xriA 1 1 e~ xr2B 
f{n 1 ,n 2 ,n 3 ,n A ,n 5 ;r,x) = / — - / — - - TZ ^- -^-^ - T ^ r - 1 - r ~^z^- (27) 

J 4vr J 4tt r 12 r 1A r 1B r 2A r 2B 

with nonnegative integers n i5 and exchange integrals which coincide with that from the previous 
section. The direct integrals are calculated as follows. When all m = the so called master 

o 

integral is given by H10H 

f d 3 r, f d 3 ri e ~ xriA p~ xr 2B 1 1 r 

f(r,x) = -^ Y 1 - L^iL (28 ) 

J 4 it J 4 77 r lA r 2B r 1B r 2A r 12 

1 
x 
and its first derivative with respect to r 

dr' F(xr')K (V2xr')-K l (^/2xr) j dr' F(xr') I (V2xr' 

(29) 



I (V2xr) dr'F(xr')K (V2xr') + K (V2xr) dr' F(xr') I Q (y/2xr') 

Jr JO 



f'(r,x) = V2 



where 

F(x) = -e x [Ei(-2 x) + Ei(-ar)] + e~ x [Ei(x) + 2 Ei(-x) - 7 - ln(2 x)]. (30) 

This one dimensional integration can easily be obtained with the sufficient accuracy and for ex- 
ample at r = 6, x = 2 the value of / is 

/(6, 2) = 0.001 790 194 708 681 916 168 871 495 878 339 876 (31) 

All integrals with higher powers of electron distances can be obtained from recursion relations, 
which were derived in Ref. II 1 ON . Since 

/(m, n 2 , n 3 , 714, n 5 ; r, x) = f(n u n 2l n 3 , th, n 5 , 1, r x) r 2+n 1+ n 2+ n3 + n 4+ n 5 (32) 

we can present formulae at r = 1. Let us introduce auxiliary functions 

y = e~ x (33) 

y x = [Ei(x) - ln(2 x) - 7] e _a! - [Ei(-2 x) - Ei(-x)] e x (34) 

y 2 = [Ei(x) - ln(2 x) - 7 ]e~ x + [Ei(-2 x) - Ei(-x)] e x (35) 

y 3 = [2 Ei(-x) + Ei(x) - ln(2 x) - 7] e' x + [Ei(-2 x) + Ei(-x)] e x (36) 

y A = [2 Ei(-x) + Ei(x) - ln(2 x) - 7] e~ x - [Ei(-2 x) + Ei(-x)] e* (37) 



then all integrals / 

/(O, 0,0, 0,0 
/(O, 0,1, 0,0 
/(0, 1,0, 0,0 

/(i,o,o,o,o 

/(O, 0,0, 2,0 
7(0, 0,1, 1,0 
/(O, 0,2, 0,0 
7(0,2,0,0,0 
7(2,0,0,0,0 

7(0,1,0,0,1 

7(0,1,0,1,0 
7(0,1,1,0,0 

7(1,0,1,0,0 
7(1,1,0,0,0 



fl-i, n 2 , n 3 , n 4 , n 5 ; x) with J2i n i <% are the following 
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(38) 


x 2x 




(39) 


2/2 

2x 3 




(40) 


(i - yo) 2 

X 4 




(41) 


2x 4 2x 3 T V^ 2 2y J x 2 




(42) 


y 2 o-i 2/4 yo (I l\ . 37' 
2x 4 2x 3 2x 3 Vx 2 2y J 2x 2 




(43) 


</o 2 -l + </o /2 + l\ 7' 
2x 4 2x 3 V^ 2 2y 7 x 2 




(44) 


Z/o,2/4 / /' (I 1\ 
x 4 2x 3 x 2 2x 2 V2x 3 x 4 ; y ° 




(45) 


_y i _%o_^ L + ^ + 27 + p + ^\ 2 _ 
x 5 x 5 2x 4 x 4 x 2 \x 5 x 4 / 


Vx 4 + 2x 2 y 


1 f (46) 


1 /l 1 \ 

x 4 Vx 4 + 2x 3 J Vo 




(47) 


2/2 2/1 2/o 

2x 4 2x 3 2x 3 




(48) 


2/2 2/o 
x 4 2x 3 




(49) 


2 M l\ / 2 l\ 2 




(50) 


l-2/o 

..4 




(51) 



X' J 



Integrals with higher powers rii are of analogous form, they are all linear combinations of f, f, 
IJq and y^ with coefficients being polynomials of 1/x. We have generated a table of integrals with 
Y^i Tii < 37, which corresponds to a maximum value of Vl = 16. It is less than VL = 20 in the case 
of James-Coolidge functions, but anyway, requires the use of octuple precision arithmetics due to 
near linear dependence of these basis functions. 

IV. NUMERICAL RESULTS 



The matrix elements of the nonrelativistic Hamiltonian between exponential functions are ob- 
tained as described by Kotos and Roothan in [16]. The resulting expression is a linear combination 
of various Slater integrals which are calculated using analytic formulas as presented in the previ- 
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ous sections. This evaluation is fast and accurate, thus allowing the use of a large number of basis 
functions. 

Eigenvalues of the Hamiltonian matrix are obtained by inverse iteration method for various 



length of the basis set. Following Sims and Hagstrom I115h we use double basis set (fl, f2 — 2) 
with two different nonlinear parameters, which were obtained by minimization of energy at fi = 
12. For the calculations with basis functions up to fi = 20, the second nonlinear parameter is 
additionally multiplied by 1 .5 in order to improve the numerical stability for the large values of fi. 
Numerical calculations with James-Coolidge functions are performed in general using quadruple 
precision arithmetics, and for checking the numerical accuracy, the one point at r = 1.4011 au is 
calculated using the octuple precision. We observe a significant loss of digits for large values of Vl. 
The quadruple precision arithmetics was not always sufficient for the whole range of internuclear 
distances. In many cases values with tt = 20, and sometimes even £1 = 19 had to be disregarded 
due to numerical instabilities in the inverse iteration procedure. Numerical results at r = 1.4011 
and r = 6 au for various sizes of basis length are presented in Table HI The most accurate result 
obtained at the equilibrium distance r = 1.4011 au is compared in Table HI] to all the previous 
results obtained so far in the literature. 

In performing extrapolation, we observe the exponential e^ 139, convergence. In other words 
the log of differences in energies for subsequent values of tt fits well to the linear function. This 
allows for a simple and reliable extrapolation to infinity. For example at the equilibrium distance 
the parameter (3 is about 0.9, and convergence in the James-Coolidge basis preserve its exponential 
behavior for all the distances. 

The calculations for the internuclear distances r = 6 — 20 au. are performed in ECA basis 
set using octuple precision arithmetics up to fi = 16. This basis is probably the most effec- 
tive one at large distances. We also observe exponential convergence for BO and even exchange 
energies, what makes quite simple the extrapolation to infinity, see Table [III] with detailed re- 
sults for various length of the basis set. For r < 12.0 results obtained with James-Coolidge 
basis are more accurate, while for r > 12.0 the ECA basis functions give more accurate en- 
ergy. Extrapolated results for the whole BO potential curve in the range 0.1 — 20 au are pre- 
sented in Table [TV] The ECA functions have the right large r behavior and thus can be used 
to study the long range tails in the interatomic interactions. Particularly interesting is the long 
range asymptotics of the exchange energy AE, which is the difference between singlet and triplet 
electronic energies. Table |Vl presents our numerical results for the exchange energy at internu- 
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TABLE I: Numerical values for BO energy obtained with James-Coolidge functions, double basis set 
(J2,fi-2) with nonlinear parameters: a% = 0.9650, a 2 = 4.6716 for r = 1.4011 andai = 0.57050, a 2 = 

2.36925 for r = 6.0 

n N r = 1.4011 r = 6.0 

-1.173 189 743 241 832 33 -0.969 640 616 224 699 7 

-1.174 345 859 790 847 31 -0.995 618 908 838 543 1 

-1.174463 692 797 208 81 -1.000 157 746 708 230 

-1.174474 90724722129 -1.000 781343 5516901 

-1.174475 84111726947 -1.000 831730125 255 5 

- 1 . 174 475 923 062 032 59 - 1 .000 835 509 296 673 5 

-1.174475 930460077 70 -1.000 835 696 691923 9 

-1.174475 931 266 449 18 -1.000 835 707 165 465 4 

-1.174 475 931 378 701 09 -1.000 835 707 596 974 8 

-1.174 475 931 396 395 07 -1.000 835 707 643 241 4 

-1.174475 931399 45343 -1.000 835 707651449 9 

-1.174475 931400035 47 -1.000 835 707653 956 

-1.174475 931 400 162 72 -1.000 835 707 654 760 4 

-1.174475 93140019706 -1.000 835 707655 025 5 

-1.174475 931400 20848 -1.000 835 707655120 8 

-1.174475 931 400 213 00 -1.000 835 707 655 156 

-1.174475 931400 215 01 -1.000 835 707655168 9 

-1.174475 931400 215 99 -1.000 835 707655175 9 

-1.174 475 931 400 216 7(3) -1.000 835 707 655 180 4(22) 
-1.174475 93140021(6) 
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clear distances r = 6 - 20 au. It was found by Herring and Flicker in [12] that the asymptotic 
limit is e~ 2r r 5//2 (7 + 0(r -1 / 2 )) with 7 = 1.636 571 460 2.... The fit to our numerical data, see 
Fig.EQ gives 7 = 1.6(1) what can be regarded as a first numerical confirmation of the Herring and 
Flicker result. When assuming their constant, the fit of first three terms in 1/ ' y/r expansion gives: 
— 1.11732/\/f + 2.13187/r + 5.169/r 3//2 , what indicates very slow convergence of this expansion 
at typical interatomic distances r ~ 10. 



TABLE II: Born-Oppenheimer potential for the H2 molecule with different types of wave functions at 
r = 1.4011 au 



Authors 



type of w.f. 



N 



energy 



1933 James and Coolidge [14] 

1960 Kolos and Roothan [16] 

1968 Kolos and Wolniewicz [17] 

1994 Rychlewski, Cencek, Komasa [18J 

1995 Wolniewicz [19] 

2006 Sims and Hagstrom [15] 

2007 Nakatsuji et al. [20] 

2008 Cencek, Szalewicz [11] 
2010 this work 
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-1.1735 


JC 


50 


-1.174 448 


KW 


100 


-1.174474983 
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700 


-1.174 475 931 


KW 


883 


-1.174 475 930 742 


JC 


7034 


-1.174 475 93139984 


ICI 


6776 


-1.174475 931400027 


ECG 


4800 


-1.174 475 931400135 


JC 


22363 


-1.174 475 931400 215 99 
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FIG. 1 : Rescaled exchange energy e 2r r 5 ' 2 AE fitted to numerical points in Table |V] in the range r 
6 — 20 au, assuming Herring and Flicker J1211 asymptotics, straight line. 
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TABLE III: Numerical values for BO and exchange energies obtained with ECA functions at r = 12 au 

n N E{ l Y,+) e 2r AE 

-1.000002 515 756 9814 720.758 075 

-1.000002 5415619799 767.006 890 

-1.000002 545 005 4934 791.784126 

-1.000002 545 635 6130 800.649 925 

-1.000002 545 838 2275 803.989 841 

-1.000002 545 922 5168 805.314 087 

-1.000002 545 955 3661 805.779 212 

-1.000002 545 965 8295 805.916 877 

-1.000002 545 968 6694 805.953 991 

-1.000002 545 969 3389 805.963 087 

-1.000002 545 9694859 805.965 246 

-1.000002 545 969 5175 805.965 766 

-1.000002 545 969 5251 805.965 907 

-1.000002 545 969 5273 805.965 954 

oo -1.000002 545 969 5279(3) 805.965 974(10) 



V. SUMMARY 

We have demonstrated the applicability of analytic formulas for two-center two-electron inte- 
grals in high precision calculations of Born-Oppenheimer potential for the ground electronic state 
of the hydrogen molecule. The symmetric James-Coolidge basis with as much as 22 000 functions 
provided energies with relative precision of about 1CT 15 for internuclear distances up to 10 au. The 
extended Heitler-London basis with about 10 000 functions give similarly accurate description at 
internuclear distances r ~ 10 au and greater. The extension of the analytic approach to excited 
states of H 2 and other diatomic molecules, such as HeH + requires the evaluation of Slater integrals 
with arbitrary nonlinear parameters (Kolos-Wolniewicz basis [13]). These integrals are solved an- 
alytically, see Ref. |i0|], but the explicit expression for higher powers of electron distances is very 
lengthy. This may need further work to find a more compact analytical form. 

Further improvement in energy requires more accurate calculations of relativistic effects, pos- 
sibly including nonadiabatic corrections. For this, integrals with quadratic inverse powers of in- 
terparticle distances have to be worked out analytically, and this is not a simple problem, at least 
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TABLE IV: Numerical values for BO potential at different internuclear distance r. Results are obtained by 



extrapolation to complete set of basis functions 



r/au 



E(r) 



r/au 



E(r) 



0.10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
1.10 
1.20 
1.25 
1.30 
1.32 
1.34 
1.36 
1.38 
1.39 
1.40 
1.41 
1.42 
1.44 
1.46 
1.48 
1.50 
1.55 
1.60 
1.70 
1.80 
1.90 
2.00 
2.10 
2.20 
2.30 
2.40 
2.50 
2.60 
2.70 
2.80 
2.90 



7.127216 731132 
2.197 803 295 226 
0.619 241659 796 
-0.120 230 341178 
-0.526 638 758 743 
-0.769 635 429 485 
-0.922 027461527 
-1.020 056 666 360 
-1.083 643 239 958 
-1.124 539 719 546 
-1.150 057 367 738 
-1.164 935 243 440 
-1.169 419 627390 
-1.172 347149 038 
-1.173138 736 333 
-1.173 734 874958 
-1.174148 498 570 
-1.174391683 632 
-1.174452917278 
-1.174475 714220 
-1.174461370 870 
-1.174411141239 
-1.174207 836 585 
-1.173 875 042 749 
-1.173 421418 292 
-1.172 855 079 578 
-1.170 994919897 
-1.168 583 373371 
-1.162 458 726 898 
-1.155 068 737611 
-1.146 850 697029 
-1.138132 957132 
-1.129163 836101 
-1.120132116849 
-1.111181765 204 
-1.102 422 606011 
-1.093 938129 955 
-1.085 791237396 
-1.078 028 484183 
-1.070 683 233 481 
-1.063 778 008 806 



0(55) 

18(33) 

226(60) 

823(15) 

001(11) 

9092(80) 

4636(26) 

5151(4) 

8343(8) 

8709(7) 

5650(3) 

3099(7) 

9022(7) 

0904(8) 

4793(8) 

3451(3) 

4193(9) 

2532(6) 

4574(4) 

4434(5) 

6800(11) 

2317(11) 

0950(11) 

2034(8) 

0817(12) 

5838(12) 

0180(6) 

4593(8) 

4588(5) 

6094(11) 

6887(28) 

6480(34) 

3193(40) 

2218(48) 

4391(17) 

326(65) 

879(75) 

1321(12) 

8287(46) 

4249(55) 

0211(65) 



3.00 

3.10 

3.20 

3.30 

3.40 

3.50 

3.60 

3.70 

3.80 

3.90 

4.00 

4.20 

4.40 

4.60 

4.80 

5.00 

5.20 

5.40 

5.60 

5.80 

6.00 

6.50 

7.00 

7.50 

8.00 

8.50 

9.00 

9.50 

10.00 

10.50 

11.00 

11.50 

12.00 

13.00 

14.00 

15.00 

16.00 

17.00 

18.00 

19.00 

20.00 



-1.057326 268 872 
-1.051333 772 268 
-1.045 799 661432 
-1.040 717365 351 
-1.036 075 395190 
-1.031858 084 855 
-1.028 046 308 379 
-1.024 618188 410 
-1.021549 795 533 
-1.018815827696 
-1.016 390 252 950 
-1.012 359 959 683 
-1.009 256 516 261 
-1.006 895 223 822 
-1.005116 006100 
-1.003 785 658 583 
-1.002 796 816 311 
-1.002 065 057 209 
-1.001525 251886 
-1.001127 880 852 
-1.000 835 707655 
-1.000 400 548 534 
-1.000197914 480 
-1.000102106147 
-1.000 055 604 973 
-1.000 032171832 
-1.000 019 781832 
-1.000 012 856 876 
-1.000 008 755 746 
-1.000 006189 995 
-1.000 004 505 989 
-1.000 003 356174 
-1.000 002 545 969 
-1.000 001529 286 
-1.000 000 960 680 
-1.000 000 625 453 
-1.000 000 419 586 
-1.000 000 288 826 
-1.000 000 203 340 
-1.000 000146 028 
-1.000 000106 740 



6617(70) 

0178(65) 

432(17) 

395(11) 

7599(41) 

0934(29) 

7348(90) 

962(10) 

649(12) 

496(12) 

6681(55) 

166(11) 

5862(34) 

7406(46) 

3838(50) 

9706(38) 

2547(33) 

7059(27) 

6137(16) 

4173(36) 

1804(23) 

5376(60) 

0381(38) 

8089(22) 

0730(4) 

8288(55) 

4911(2) 

8268(2) 

0515(1) 

1069(1) 

4362(1) 

5754(1) 

5285(1) 

6698(1) 

7911 

6319 

3122 

2392 

5059 

2368 

1283 



for arbitrary nonlinear parameters. 

Even more important is the extension of this approach to few and many electron diatomic 
systems. Integrals with (uncorrelated) Slater functions with arbitrary nonlinear parameters can 
be evaluated using recursion relations of Ref. IllOn . this includes also the exchange integrals in- 
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TABLE V: Exchange energy AE = E u — E g at different internuclear distances 

r e 2r r~ 5 / 2 AE 

6 1.886 757 524(7) 

7 1.798 767166(8) 

8 1.736 967949(9) 

9 1.692 575 090(10) 

10 1.659 770 272(12) 

11 1.634 919 801(15) 

12 1.615 710 667(20) 

13 1.600 619 338(31) 

14 1.588 607 797(56) 

15 1.578 94736(11) 

16 1.57111318(24) 

17 1.564 718 78(50) 

18 1.5594741(11) 

19 1.5551578(22) 

20 1.551599 0(44) 

volving Coulomb interaction between electrons. More challenging would be implementation of 
coupled-clusters with the explicit r\i factor. The analytic form for the corresponding three-electron 
integrals is however not yet known. 
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